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Abstract- In this paper we study the spatial and temporal behaviour of the dy- 
namic processes in porous elastic mixtures. For the spatial behaviour we use the 
time-weighted surface power function method in order to obtain a more precisely 
determination of the domain of influence and we establish spatial decay estimates of 
Saint -Venant type with time-independent decay rate for the inside of the domain of 
influence. For the asymptotic temporal behaviour we use the Cesaro means associ- 
ated with the kinetic and strain energies and establish the asymptotic equipartition 
of the total energy. A uniqueness theorem is proved for finite and infinite bodies 
and we note that it is free of any kind of a priori assumptions of the solutions at 
infinity. 



1 Introduction 

Various theories have been proposed in literature for describing the behaviour of the 
chemically reacting media (see, for example, Truesdell and Toupin |]], Kelly @, Eringen 
and Ingram || |], Green and Naghdi || ||, Miiller Q, Dunwoody and Miiller ||, Bedford 
and Drumheller |J, etc). 

*Work performed in the context of the COFIN MURST, Italy, "Mathematical Models for Materials 
Science" and under auspicies of G. N. F. M. of the Italian Research Council (C. N. R. ). 
tTel.: +39-89-964252, fax: +39-89-964191 E-mail address: passarella@diima. unisa. it 
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Recently Ie§an |TD| has developed a theory for binary mixtures of granular materials 
in Lagrangian description in which the independent constitutive variables are the dis- 
placement gradients, displacement fields, volume fractions and volume fraction gradients. 
The theory takes into account the results established previously by Nunziato and Cowin 
fTT| , Goodman and Cowin |12 and Drumheller |13|]. The intended applications for such 
a theory are to granular composites, solid explosives and geological materials. 

In [ID] a linear theory is also presented and some uniqueness results for bounded bodies 
are established for the linear dynamic theory with no definiteness assumptions on the 
elasticities and without any restriction on the initial stresses. 

The present paper studies the spatial and temporal behaviour of the solutions to the 
boundary-initial value problems in the linear dynamic theory of porous elastic mixtures 
as developed in [p~0[1 . 

For the spatial behaviour of the dynamic processes in porous elastic mixtures we use 
the time-weighted surface power method developed in fH |. Thus, we introduce a time- 
weighted surface measure associated with the dynamic process in question and then we 
establish a first-order partial differential inequality whose integration gives a good in- 
formation upon the spatial behaviour. Then we obtain a more precisely version of the 
domain of influence in the sense that for each fixed t e [0, T] the whole activity is van- 
ishing at distances to the support of the given data on [0, T] greater than ct, where c is a 
constant characteristic to the elastic mixture. A spatial decay estimate of Saint-Venant's 
type is established for describing the spatial behaviour of the dynamic process inside of 
the domain of influence. 

As regards the temporal behaviour of the dynamic processes in porous elastic mixtures, 
we introduce the Cesaro means of various energies and then establish the relations describ- 
ing the asymptotic equipartition of energy. In this aim we use some Lagrange identities 
and the method developed by Day JL5|] and Levine |16 . 



The plan of our paper is the following one. In the Section 2 we present the basic 



equations of the linear dynamic theory of porous mixtures developed in JT0|. Some con- 
stitutive assumptions and other useful results are also presented. The auxiliary identities 
are established in the Section 3, while in the Section 4 a time-weighted surface measure is 
defined and its properties are studied. Moreover, a first-order partial differential inequal- 
ity is established for this measure. The main result concerning the spatial behaviour is 
presented in Section 5 and some uniqueness results are obtained as a direct consequence. 
In the Section 6 we introduce the Cesaro means of various energies and establish the 
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asymptotic equipartition of the total energy. 



2 Basic equations 

Throughout this article, the motions of continuum are studied respect to a fixed ortho nor- 
mal frame in IR 3 . Then, we deal with functions of position and time. Moreover, it is useful 
stress that in the following text the tensor components of order p > 1 will appear with 
Latin subscripts, ranging over the integers {1,2,3}, and summation over repeated sub- 
scripts will be implied. Greek indices are understood to range over {1, ... ,9} if they are 
lower case letters, or over {1,2} if they are upper case letters; summation convention 
is not used for these indices. Occasionally, we shall use bold-face character and typical 
notations for vectors and operations upon them. Superposed dots or subscripts preceded 
by a comma will mean partial derivative with respect to the time or the corresponding 
coordinates. 

Let B be a bounded or unbounded regular region in the physical 3-dimensional space, 
whose boundary dB is a piecewise smooth surface. A chemically inert binary mixture of 
two interacting porous elastic solids, C\ and c 2 , in a given reference configuration, is into 
B. 

The positions of particles of c\ and c 2 at time t are x and y respectively, i. e. 

x = x(X,£), y = y(Y,t) X,YeB, tel, 

in which X and Y are reference positions of these particles, I = [0, oo). By following 
Bedford and Stern [T^J, we assume that X = Y. 

Let the top label a refer the various fields to the constituent c a . Taking into account 
the linear theory, the behaviour of a binary mixture of elastic solids is governed by the 
local balance equations (see Iesan jlQ|| ) 

(i) 

h\f + g (a) + Q (a) t {a) = P {a) X (a) (p {a \ on B x (0, +oo) 

In these equations, S^ a \ are the stress tensor and body force associated to c a ; p is 
the vector field for characterising the mechanical interaction between the constituents c\ 
and the equilibrated stress vector, intrinsic and extrinsic equilibrated 

body force associated to c a , respectively. 
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Moreover, is the displacement vector fields associated to c a \ tp^ is the changes in 
volume fraction starting from the reference configuration to c a . 

Finally, , X are the bulk mass density and equilibrated inertia of the material c a 
in the reference state. 

According to classical interpretation of system (|IJ), we assume that 

i. v!f\^ eC 2 > 2 {BxI); 

ii. S£\h^ E C^{B x I), pi £ C°>°(B x 7); 

iii. ft\g {a \£ ia) e C°'°(B x /), P ( a \x (a) e C°{B), 

where, B is the closure of B. 

Then, we introduce the 29-dimensional vector field 



E(U) = {e u -(U),^(U),^ 1 ),^,rf l (U),^ ) ,^ ) }, 



with 



and 



eij (U) = \ («g + uj J ) , 9ij (U) = u# + u# , 



W . .,(2) 



on B x I. 



di(U) = 4 1} - uf , 
Now, we define the magnitude of E(U) by 

|E(U)| = {E h(U)ey(U) + ^(U)^-(U) + <^(U)^(U)+ 



(2) 



a=l 



1/2 



+4(U) + ^ ) (U)^ ) (U)]} 

Our attention is focused on homogeneous, centrosymmetric mixture, by supposing the 
initial continuum is free from stresses. 

In the context of our theory, the internal energy density associated to U is given by 
1 



W(V) = - [A,>^(U)e rs (U) + C ijrs9ij {U) 9rs (U) + C^ (1 V (1) + 
+/W 2 ) + atijtpftpf + Ivvfvf + a^(UK(U) 
+B ijrs e ij (V)g rs (U) + D tj e tj (U)(p^ + E ij e ij (\J)^ + 



+ 



(3) 
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The material coefficients, appearing in previous equations (^]), are constants and they 
obey the following symmetry relations: 



(4) 



(5) 



■A-ijrs -^-jirs -^-rsij ; £$ijrs -^jirs 5 ^ijrs ^rsij ? ®ij dji 5 
OCij lij Ijii ^jii ^ij ^ji- 

The constitutive equations are 

Sji (U) = (Aji rs + B rs ji)e rs (XJ) + (Bij rs + Cji rs )g rs (U) + 

4 2) (U) = B rsij e rs (U) + C ijrs g rs (U) + M ijV ^ + N^ 2 \ 
g (i)(lJ) = -D rs e rs - M rs g rs - Cv {1) - V 2 >, 
0<2)(U) = -E rs e rs - N rs g rs - ry>« - 
Pi(XJ) = ciijdjiXJ) + bijtpy + Cijip ( f , 
h^(U) = auwf + /3 ij( pf + bpdjiV), 

hf\u) = Pjwf + 7ij <pf + c^U). 
Let A x = \\a KL \\ (K,L = 1,...,20) 

arA = ^4rA, a r (9+A) = B T a, Sri9 = D r , a r2 Q = E F , 

fl(9+r)A = B Ar , a( 9+ r)(9+A) = Cta, «(9+r)i9 = M r , a( 9+r )2o — N r , 

dig A = -Da, Ol9(9+A) = Ma, a.1919 = C; O1920 = t, 

&20A = Ea, O20(9+A) = ^Va; 5-2019 = T i 5-20 20 = A 4 ) 

where we have called the nine index combinations (ij) or (r s) by capital greak letters 
(z.e.r , A , and so on). Now, let O be the empty matrix 20 x 9 and A2 = ||&kx|| (K, L = 



(6) 
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1,...,9) be 



i(3+i) u iji 



b(3+i)j - bji, &(3+i) (3+j) 



a 



i j ■ 



A(6+j) c iji 



b(3+i) (6+j) — Pij, 



(7) 



b(6+i)j — Cji, &(6+i) (3+i) — (6+j) — 7ij- 

Then, the energy density (|3]) assumes the form 

1 29 - 1 
W(V) = - ]T A K l^(U)£7 l (U) = -E(U).^E(U), 

1 K,L=1 Z 

where the matrix A = \\Akl\\ (K, L — 1, . . . , 29) is defined by 

A! O 



A 



O 1 



A 2 



(9) 



and O t is the transposed matrix of O. 

In what follows, we assume that p^ a \ X are strictly positive and VF(U) is a positive 
definite quadratic form; thus, there exist the positive constants £ m and £m so that 

UE(U)| 2 > 2W(U) > £m|E(U)| 2 (10) 

where £ m is the minimum elastic modulus and £m is the maximum elastic moduli. 
Let 

S(U) = {Sff (U) , 4 2) (U) , (U) , g® (U) , Pi (U) , (U) , hf ] (U) } , 
then the magnitude of S(U) is defined by 

|S(U)| = {E[sjl?(V)Sj?(V) + Ji a \U)hP(V) +g^(V)g^(V)} + ft (U)p,(U)} 1 / 2 

a=l 

Taking into account the equations (f| H, 0; & ED) we prove that 

|S(U)| 2 = AE ■ AE = E • ^ 2 E < £ M E ■ „4E = 2£ M W(U). (11) 
The surface tractions s(°)(U) and h^ a \\J) are defined by 



s\ a \v) = 4 Q) (U)n i; fr<«>(U) = /4 Q) (U)n„ 



(«), 



(12) 



where n is the outward unit normal vector to boundary surface. The relations (|TT|, [T2|) 
imply that 



^[4 a) (U)4 a) (U) + ^(U)^(U)] < |S(U)| 2 < 2CmW(U). 



a=l 



If we introduce the notations 



A.. - C- 



a . . — E- ■ 4- N- ■ 



Tij = Dij + My 

the equations (|3p become 

W(U) = \ [a^ugllg + dyr.ugug + C^ (1 V (1) + /W + 

+a ii ^Vj ) + 7ij^Vf + 0^^(11)^(11) + b ijrs uf] uJ5+ 
+r,, W £V (1) + ^'4^ (2) + ^ugV W + ^£V 2) + 

Using the symmetry relation (|]) ; we get 



Qijrs &rsij j 



(13) 



(14) 



(15) 



(16) 
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The constitutive equations (0) become 

Sff(U) = a ijrs u^} + b ijrs uf) s + + a ijV {2 \ 

4 2) (U) = b rsijU W + d l3TS uf} + MijipW + N^ 2 \ 

<?«(U) = -r rsM « - M rs uf} - - V 2) , 

5 (2) (U) = _ CTrsU (l) _ ^(2) _ r ^(l) _ ^( 2)) 



(17) 



Pi(U) = dijdjiXJ) + b^cpy + Cij(p[ 



,( 2 ) 
3 ' 



/if >(U) = Oi^y + Pijcp)? + bjidjiU), 



(2) 



hf \v) = Pjwy + jiw)- 1 + ^(1]). 

It follows from the equations (0), (0), (0) that 



(2) 



2W(V) = E[4 a) ( U )4? +Pi(U)^(U) + ^ a) (U)^ } -0<«>(U> 



(a.) 



a=l 



and 



w(u) = £ [4 a) (u)«g } + Pi (u)d i (u) + ^ a) (u)^ - ^ Q )(u)0 



(If 



(19) 



We consider the initial-boundary value problem V defined by the equations of motion 
(0), the geometrical equations (Q) and the constitutive equations @ and the following 
initial-boundary conditions 



(a) (a) . (a) . (a) 



(20) 



onBx {0}. 



and 



(a) ~(ol) v=i t 

u\ —u\' on Si x I, 



(a) ~(a) T 

s\ — s\ ' on S 2 x I, 



pi") = (p( a ) on S 3 x /, = on S 4 x I, 

where Sj (? = 1, 4) are the subsets of 55 such that 



(21) 



Si u s 2 = s 3 u s 4 = 95, Si n s 2 = s 3 n s 4 
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The terms on right-hand in equations ( |2D| ) and (ET) are prescribed continuous functions; 
along with f M , f ( 2 ) , i^ l > , £( 2 ) these constitute the external data of the problem V. 

An array field U = {u«, \P\ tpM, c^ 2 )}, meeting all equation (g, |, |, H) and (P), 
will be referres to as a (regular) solution of the problem V. 



3 Auxiliary identities 

In this section we establish some integral identities that we will use in next sections. 

Lemma 3.1 Let U be a solution of initial-boundary- value problem V. Then, for every 
regular region P C B with regular boundary dP, it follows that 

- / e- xt y[p ia) u^\t)4 a \t) + p (Q) X (a) ^ a \t)^ a \t) + 2W(U(t))]dv+ 
2Jp *=i 



+ | f I e~^j:[p^u { t\s)v!t\s) + P (a) X ia) ^(s)^(s)+ 
2 Jo Jp ^ 



+2W(V(s))]dvds = T / e" As j2lP (a) ft\s)u^\s) + p {a k ia \s)<p {a \s)}dvds+ 
Jo Jp a=1 

+\ I E[/° (a) «i a) ( )^ Q) ( )+P (a) * (a) ^ {a \0) 2 + 2W(V(0))]dv+ 

^ J P 1 

f I e- Xs J2[sP(V(s))u^(s) + h^(U(s))^ a \s)]dads, 
Jo JdP ^ 



ft 

+ 



(22) 

where A is a positive parameter and t £ I. 
Proof. The equations (|l|) and (|T9"D lead to 

d rl 2 



a=l 
2 



e~ Xs 



a=l 

+ [4 a) (U( S ))^^) + hf(V(s))^\s)]. . 

By an integration of the equations (^) over P x [0, t] and by using the divergence theorem, 
we obtain the equation (P2"|). ■ 
If we introduce 

E(t) = I \ Y.\P {a)ibi t\t)u ( t\t) + P (a) X {a) ^ a \t)^ a \t) + 2W(V(t))]dv, (24) 

JB 2 a=l 
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then, for A = and P = B the equations (E2p reduce to 

ft r 2 



S(t) = £{0) + f f Y[p (a) ti a \s)u\ a) {s)+p {a k {a \s)v {a) {s)]dvds+ 

Jo JB a=1 



+ 



f I EK (a) (U(s))«! Q) ( S ) + h^(U(s))^(s)]dads. 

JO JdB^ 



(25) 



We note that S(t) is a measure of the energy stored in B at time t. 

Lemma 3.2 Let U be a solution of initial-boundary- value problem V. Then, for every 
regular region P C B with regular boundary dP, it follows that 



a=l 

= T / E{[P (Q) ^ a) ( s )^ (a) (s) + P (a) X {a) ^ a \s)^ a \s)]-2W(U(s))}dvds+ 

+ / Eb (a M a) (o)wS Q) (o) + p (q) x (Q) ^ (Q) (o)^ (Q) (o)]^+ 

+ f I E[4 Q) (U( S ))«! a) («) + h^{\J{s))^\s)]dads+ 

JO JdP a=l 

Proof. The relations ([!]) and (|18D imply that 

^[jl[p {a) ut\s)ut\s) + pM X {a) ^\s)^(s)} = 

as a =l 

= Y,[p {a) u ( r\s)u ( r\s)+p {a) x (a) ^\s)^\s)]+ 



(26) 



t e I. 



a=l 
2 



(27) 



+ E [ P {a) ft\s)uf(s) + p^^\s)^ a \s) - 2W(V(s))+ 

a=l 

+[4 a) (u( S )) W s a) ( S ) + / i f(u( S ))^( S )],; . 

The relation ( p6| ) follows from ( p7|) by an integration of over P x [0,t] and by using the 
divergence theorem. ■ 

Lemma 3.3 Let U be a solution of initial-boundary-value problem. Then, for every 



10 



regular region P C B with regular boundary dP, it follows that 

2 



2 / Y.[p {a) ^\t)ut\t)+P {a) X (a) <p<«\t)<pM(t)]dv 



a=l 

2 



= I Eb (a) ^ ) (o)«! Q) (2t) + P (a) ^ ) (o) M ! Q) (2t)+ 

+p (a) X {a) v? (a) (0)^ {Q) (2t) + p {a) (^ (a) (2t)^ (Q) (0)]^+ 



+p {a H i - a \t - s)v {a \t + s) - p {a) &\t + s)^ a) {t - s)]dvds+ 

+ tf EK (Q) (U(t - s))u^\t + s)- s\ a) m + s))uf(t - s) + 

JO JdP^ x 

+/i (a) (U(t - s))^ {a \t + s) - h {a) {U{t + s)) V {a \t - s)]dads, t e 

Proof. For every function <fi G C 2 (J), it is holds the following identity 
<j>(t - s)(p(t + s)- '(j){t + S)(j){t - s) = 

-t{# - s )^(* + ") + # - *)# + ")}« se[o,t], tei, 

as 

In view of the relation ([I]), we have 

E [P {a) ^\t - s)uf(t + s)- p^ut\t + s)v!f(t - s)] = 

01=1 

j:{p {a) fl a \t - S)u^(t + 8)- P ^ft\t + S)vi a \t - S) + 



a=l 

+[s£>m - s)) U r(t + S )- spm + s)) U r{t - 

+ [S£\U(t + 8))u^{t -8)- - 8 ))v${t + S)]} + 

- Pi (U(t - s))di(U(* + s)) + Pi (U(t + s))d,-(U(t - s)), 
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and 



i2lP [a) X H ^ {a \t-s)ip^\t + s) -p {a) X {a) 0W(t + s)<pM{t-s)] 



a=l 

2 



Y^{p {a) £ {a) (t - s)<p( a \t + s)- p {a) £ {a) (t + s)v {a) (t - s)+ 



a=l 



+[h<f\v(t - s))ip {a \t + s)- hf(V(t + s))<p<- a \t - s)]j+ (31) 



+[h) a) (U(t + s))<p$>(t -a)- h]>(XJ(t - s))<p%'(t + s)]+ 

+g(U(t - s))<pM(t + s)) - g(U(t + s))<pW{t - s))\ 
Further, with the help of (|Hi|)-(|i7|) we prove that 

Y.{sfm + s))v${t -a)- S$?m - 8 ))v$(t + a)+ 

a=l 

+[hf{V{t - s))<p<- a \t + s)- hf(U(t + s))<pM(t - s)] tJ + 

+[h<f l) (U(t + s)) V f(t -s)- hf{XJ(t - s))<pf(t + s))+ ( 32 ) 

- Pi (U(t - s))di(U(t + s)) + Pi (U(t + s))di(V(t - s))+ 

+g(U(t - s))<pW(t + s) - g(U(t + s))<p&>(t - s)} = 0. 
Then, the equations (^9|)-(p^) imply that 



EiP^f^it - S^it + 8)- P ^fP(t + 8)li a \t -S) + 



a=l 



(33) 



p {a k [a \t - s )ip^(t + s)- p^k {a \t + s)v [a) (t - S)+ 

+ [hf\\J(t - s))^\t + a) - + 8))<p<*\t - 8)]j} = 

= -j- i:{p {a) ^\t - a)v<r\t + S )+ p^ufit - s^it + S ) 

as a=l 

The equation fl28|) is reached by performing an integration of the equations (^) over 
P x [0,t] and then by using the divergence theorem. ■ 
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4 A time weighted surface measure 



Having fixed a time Tel, for the external given data of the problem V we define the set 
D T by: 

i. if xG B, then 

a! 1} (x) ^ or of } (x) ^ or af } (x) ^ or df } (x) ^ 

or 

rf } (x) ^ or v4 1} (x) ^ or ^ 2) (x) ^ or ^ 2) (x) ^ 
or there exists r G [0, T] such that 

/f(x,T) ^ or /f(x,r) ^ or £«(x,r) ^ or ^(x,r) ^ 0; 

ii. if x G dB, then there exists r G [0, T] such that 

4 1} (x,rK (1) (x,r)^0 or Si (2) (x, r)nf»(x, r) ^ 

or 

/i (1) (x,r)v? (1) (x,r) ^ or /i (2) (x, r)y? (2) (x, r) ^ 0. 

The set .Dy represents the support of the initial and boundary data and the body force 
on the time interval [0, T]. In what follows, we assume D T is a bounded set. 
We consider a nonempty set D T which is such that D T C D T C B and 

i. if Dt D B ®, then we choose D T to be the smallest bounded regular region in B 
that includes D T ; in particular, we set D T = D T if D T it also happens to be a regular 
region; 

ii. if 7^ D T C then we choose D T to be the smallest regular subsurface of dB that 
includes D T ; in particular, we set D T = D T if D T is a regular subsurface of dB; 

iii. if Dt = 0, then we choose D T to be an arbitrary nonempty regular subsurface of 
dB. 

Now, we mean the set D r , by 

D r = {xGS : D;nS(x,r) ^ 0}, r > 0, (34) 

where S(x, r) is the closed ball with radius r and center at x. Further, we use the notation 
B r for the part of B contained in B \ D r and B(ri, r 2 ) = B r , 2 \ B ri , r\ > r 2 ; S r denotes 
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the subsurface of dB r contained into inside of B and whose outward unit normal vector 
n is forwarded to the exterior of D r . Of course, taking into account that for each r > 0, 
D T C D r and D T n B r = 0, we get 



ftW = 0, v<l a) = 0, = 0, C (a) = on B r , 

f^ = 0, £(<*) = onB r x[0,T], (35) 

j a >4 a) = /i( Q )0 (a) =0 on (B r n SB) x [0, T]. 



For a fixed positive parameter A and for any r > 0, t G [0,T], we define the time- 
weighted surface power function P(r,t) 

P(r,t) = -[*[ e- Xs E[ fl< (a) (U( fl ))ui a) ( fl ) + hM{V( 8 ))<pV>{8)] dads. (36) 

JO JS r ^ 

In the following Lemmas, we show some relevant properties of the function P(r,t). 

Lemma 4.1 Let U be a solution of initial-boundary- value problem V and be the 
bounded support of the external data on the time interval [0, T]. Then, the corresponding 
time-weighted surface power function P(r, t) is a continuous differentiable function on 
r > 0, t G [0, T] and 

^P{r,t) = -f e- xt j:[sl a \v(t))u^\t) + h^(V(t))^(t)}da-, (37) 

OZ JS r a=l 

y r P(r,t) = -\ \ s e- xt E [p (a) «i a) (*)«i a) (*) +P (Q) * W (Q) W^ (Q) (t) + 

+2W(U(f))l da-\C \ e-" s £ (38) 

Z JO J-SV- a=l 

+p (a) v? (a) (s)v9 (a) (s) +2W(U(s))] dads. 

Moreover, at a fixed t G [0, T], P(r,t) is a non-increasing function with respect to r. 
Proof. The equation (^) is an immediate consequence of the definition of P(r,t). 
The Lemma 3.1 for B(r\,r2) and the divergence theorem imply 
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P(n,t) - P(r 2 ,t) 



t r 2 

e~ Xs 



10 JB(n,r 2 ) a=1 



Eb (a) /i a) (^f'(s) + P (a) £ {a) (s)^ a \s)} dvds+ 

(39) 

2 Jo JB( ri ,r 2 ) OS 

+2W r (U)(s))] duds, 

with < r2 < ri and fixed £ > 0. By using equations (|35"D and performing an integration 
by parts, the relation (RSf) becomes 



P(n,t) -P(r 2 ,t) = ~f e~"j:\plM a \t)<£ a \t) +P ia) * {a) ^ a \t)^\t) + 

1 JB( ri ,r 2 ) a=l 

+2W(U(t))]dv-^ f f e- As ^[pM«; Q) ( S )^ (Q) ( S ) + p^ X ia) ^ a \s)^ a \s) + 

2 Jo JB(ri,r 2 ) a=1 



+2W{\]{s))]dvds. 

(40) 

This relation straighty leads to (38). 



Since we have assumed that p( a \ are strictly positive, A is positive and W(U) is a 
positive definite quadratic form, the equation (|40|) gives 



P(n,t) < P(r 2 ,t) with n > r 2 . - (41) 



Lemma 4.2 Let U be a solution of initial-boundary-value problem V and D? be the 
bounded support of the external data on the time interval [0,T]. Then, the function 
P(r,t) satisfies the following first-order differential inequalities, for any r > 0, £ G [0,T] 



where 




b.\p( rt t)\ + JLp(r t t)<0, 
c or 



0_ 

dt 



P{r,t) 



+ ^P(r,t)<0, 



^- W ith m = min {p« p( 2 >, p« x « p ( 2) x ( 2 )}. 



(42) 
(43) 

(44) 



and £m is the maximum elastic moduli. 
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Proof. It follows from Schwarz's inequality and the arithmetic-geometric mean 
inequality 



E[4 Q) (u(t))^ ) (t) + / i (a) (u(t))0 (Q) (t)" 

a=l 



(45) 



a=l 



where e is an arbitrary positive constant. 



sp^ X 



(a) 



/i (a) (U(t))/i (a) (U(t)) 



Using the relations (p6| , 45] , [|4 , [13| ) and taking e = c, we deduce 



Z JO J S r a=l 



01=1 

+2W(U(s))]dads 



(46) 



for r > and < t < T. 

Similarly, from (j37[) we obtain 



P{r,t) 



< - I e 

~ 2Js, 



-xt 



a=l 



(47) 

+2W(U(t))]da. 

for r > and t € [0, T}. By the equation (|38|) and the relations (|46|) and (f47|) we obtain 
the result. ■ 

Lemma 4.3 Let U be a solution of initial-boundary-value problem V and Dt be the 
bounded support of the external data on the time interval [0, T]. Then, it follows that 



moreover 



P(r, t) > 0, for r > 0, < t < T : 
P(r,t) = E{r,t) , 



(48) 
(49) 



where 



P(r,t) = ~ / e- xt £ [p (a) u^(t)u^(t) + p^ X {a) ^ a \t)^ a \t) + 

2 J B r a=1 



+2W(V(t))} dv + ± ft e-" s <T [ p W«S°)( a )u i W( fl )+ 



(50) 
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Proof. If B is a bounded body, then the variable r ranges on [0, L], where 

L = max {min{ [(xi — Ui)(xi — y^)] 5 : y G D^} : x G B } < oo. (51) 
Starting from the definition of Dj< and by using the relation fl36|), we obtain 

P(L,t) = 0, 0<t<T; (52) 



thus, the equation (|4l| ) implies the relation fl48|). 

If B is an unbounded body, then the variable r ranges on [0, oo). The inequality (|43| ) 
is equivalent to 

^! F(r ' t)+ l FM) -°' (53) 



and 



l§i P ^ t)+ W/^-^ (54) 



If we choose the initial condition (ro,to) such that to G [0,T] and tq > eto and we put 

t = to H in inequality fl53f ), then 

c 

^[P(r,t + ^)]<0 ; (55) 

thus, 

P(r, t + ^— ^) < P(n, t + ri ~ r ° ) with r > r x . (56) 
c c 

For r = r and r\ = r — cto, we get 

P(r ,to)<P(ro-ct ,0). (57) 
Similarly, by setting t = t in fl54|), it follows 



±[P(r,t - r —^)}< 0, (58) 
ar c 

so that 

P(r + ct ,0) < P(r ,t ). (59) 

Taking into account P(r — ct ,0) = and P(r + ct o ,0) = 0, the relations (|5TD, (|59|) 
imply that 

P(r ,t ) = 0. 
17 



Of course, for r — > oo in the above relations, it follows 



Pfoo,< 



lim P(r ,* ) 

r — >oo 



0. 



(60) 



and, by (f|3), we conclude that the relation fliq) is true. 

The equation (|4*9|) follows from the relation (EDI) by means of the use of the relations 
52]) and (BO 



5 Spatial behaviour 

By the properties of the time-wieghted surface power function P, we establish the theorem 
that gives a complete description of the spatial behaviour of the elastic process in question 
outside of the support of the external data. 

Theorem 5.1 Let U be a solution of initial-boundary- value problem, Dt be the 
bounded support of the external data on the time interval [0, T] and let P(r,t) be the 
time-weighted surface power measure associated with U. 

i. Spatial behaviour: For each fixed t £ [0, T] and < r < ct, we have 

P(r,t) < P(0,t)exp(--r); (61) 

c 

ii. Domain of influence results: For each fixed t G [0, T] and r > ct, we have 



(/ ,(i) = 0, ipW = onB r x[0,T]. 



Proof. The equations (|42|), fl4|) give 



so that we obtain the equation (pT|) . 



(62) 



d X 

— \exp(-r)P(r,t)}<0, 0<r<ct, 0<t<T, (63) 
or c 



If we choose t £ [0, T] and we set r = ct in (^), then 

|[P(cM)]<0, (64) 

and 

P(r', t) < P{ct, t) < P(0, 0) = 0, V r' > ct. (65) 
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In order to (fnf), we have for all r > ct 

P(r,t)<0, (66) 
and, taking into account (|48|) , we obtain 

P(r,t)=0. (67) 
Now, the equations (|6"TD, (f49|), (|50|) imply 

2 j-Br rri L 



+2wm))} dv+^-ff e - xs x; [p (q) ^ ) (^)^ ) (,)+ (K 



+p (a) ipW(s)tpW(s) + 2W(U(s))\ dvds = 0. 



Since and X are strictly positive, A is positive and W(U) is a positive definite 
quadratic form, we have 

1^ = 0, mS 2) = 0, (^ (1) = 0, vi (2) = on B r x [0,T]; (69) 

so, by (|35|) , we obtain the equations (|62]) . ■ 

If we put t = T and r = cT in the equations (|62|) for t = T and r = cT they imply that 
the set -D C T covers a domain of elastic disturbances produced by the the data at time T, 
i.e. 

vjp = 0, uf ] = y? {1) = 0, <p {2) = on S cT x [0, T}. (70) 



This result is known as a so-called domain of influence theorem (see Gurtin ||18|| ) 

As an immediata consequence of the equation fl6"2|), we establish the following uniqueness 
result valid for a bounded or unbounded body 

Theorem 5.2(Uniqueness) It exists at most one (regular) solution for the boundary- 
initial- value problem. 

Proof. Thanks to the linearity of the problem, we have only to show that the null 
data imply null solution. Let U = {u*- 1 - 1 , u^, (p^\ fi^} a solution corresponding to null 
data. Since the set = for each T G (0, +oo) and the function P(r,t) = 0, we can 
conclude that 

uf ] = 0, ui ] = 0, = 0, <p {2) = on B x /. . 
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We remark that if B is a bounded regular region, for values of T big enough, then it 
exists a value of £ £ [0, 2^] having the property that D ct = B, the relation fl62|) becomes 
superfluous and the behaviour of solutions is full described by the relation (|61~D. On the 
other hand, for values of T sufficiently small, the behaviour of solutions will be described 
by the relation (|62|) almost as in B. Similar arguments are valid for an unbounded regular 
region. 



6 Asymptotic equipartition of energy 

Throughout this section we study the time asympototic behaviour of the solutions of 
the initial-boundary value problem Vo for the bounded regular region B defined by the 
following equations of motion 

Sjg + (-l)°p f = pW«i a) , 

(71) 

h\f + g ( a ) = p( Q ) x {a) (p {a \ on B x (0, +oo), 

the geometrical equations (0) and the constitutive equations (||), the initial conditions 
(|20| ) and the boundary conditions 

uf ] =0 on Ex x /, s[ a) = on S 2 x I, 

(72) 

= 0, on S 3 x /, = on S 4 x /. 

Now, we introduce the Cesaro means of various energies associated with the solution 
U of the problem Vo : 

KUt) = E tJt f I P {a) X (a) <pW{ 8 )<pW(8)dvda, (73) 
£r 1 It Jo Jb 

S c(t) = E 7 T / W(V(s))dvds, 



and 

Kc(t)=K%(t)+K%(t), (74) 



20 



If measY,! = , then there exists a family of rigid motions and null change in volume 
fraction which satisfy the equations (|7T| , 0, [5], |20[) and (^). We decompose the initial 
data af^ and as 

aS a) = aS Q) +4*\ dl a) = ^ (a) + 4 a) , (75) 

where af^ and a/ are the rigid displacements determined so that A[ a ^ and A^ satisfy 
the normalization restrictions 

/ p {a) A i t ) dv = Q, f p {a) e ljk x j A i k a) dv = 0, 

J B J B 

p^A\ a) dv = 0, I p {a) s ljk x j A k a) dv = 0, 



(76) 



IB JB 

and Eijk is the alternating symbol 
We put 



and 



and 



C\B) = {v whit Vi E C l (B) : v { = on E : if meas^ ^ 0, 

f 2 f 2 

or P^Vidv = , 1^2 P^' ] e ijkXjVkdv = 0} if measEi = 0, 

C 1 (B) = {CeC 1 (B): C = 0onS 3 }, 
Wx(_B) = the completion of C 1 (S) by means of || ■ ||wi(s), 



and 

W\(B) = the completion of C X {B) by means of || ■ ||wi(.b)- 

The spaces W m (B) represents the familiar Sobolev space and W m (B) = [W m (B)} 3 . 
The equation (|T0|) assures that the following Korn's inequality [|19| holds 

2 2 

f ^2W(V)dv> mi f J2(vl a) vl a) +X (a) (p (a) (f) ia) )dv, mi = const. > 0, (77) 

JB a=l JB q=1 

for every V = { v« , v^ 2 ) , <jF> , ^ } : EW X (B), ^ G W x (B) . 

If measLi = , then we shall find it a convenient practice to decompose the solution 
U = {u^ , vS 2 \ <£>M, y}( 2 )} of the problem Vq in the form 

= df + thf + vV , ^ Q) = (Q) , (78) 
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where V == {v^, v( 2 \ <p^} G W X (B) x Wi(B) x Wi(.B) x W^B) represents the 
solution of the problem Vo with the initial data {af \ and {af \ fi^}. 

Theorem 6.1 Let U is the solution to the problem Vq. Then, for all choise of initial 
data with at"' G Wi(B) , a( Q ) G W (B), G , ^ G W (5). Then the 

following asymptotic behaviour of the solution U holds: 

i. if measSx 7^ , we have 

lim /C c (t) = lim S c (t) = ls(0); (79) 

ii. if measEi = 0, we have 

lim Ka(t) = lim S c (t) + ~ f £ P (a) af dv = 
t~*oo t^oo z Jb ~, 

(80) 

= \m + \ i E P^^a^dv, 
2 2 ^ 

where is defined by fl2"4|) . 

Proof. Respect to the problem P it follows that fW = 0, f (2) = 0, = 0, £ (2) = 
and the boundary conditions (f72|) are verified. Thus, the equation becomes 



£(t)=£(0), t>0. (81) 

so that 

JC c (t) + S c (t) = £(0), for all t > 0. (82) 
On the other hand, the relations (p6|) and (p8|) imply 



K c {t) - S c {t) = -1 jf Ef2p (a) ^ (Q) (0)^ (Q) (0) + 2pW ^(0)^(0)+ 



Q = l 



+p( Q) {^ (a) (0)^ (a) (2t) + ii^ (0)wS a) (2t) + X W <^(0)^(2i)+ (83) 
+p (Q) ^ (a) (2t)y? (cl) (0) W t > 0. 



The relations flT0|), (g§ and (|T|) imply 

/oW«S a) (s)^ o) (s)dt; < 2£(0), 



p (a) x («) (pW^tpW^dv < 25(0), (g 4 ) 

/ ^(s)^(s)^ < ^ / E W(U(s))^ < ^-£(0), 
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thus 

lim - / p^u\ a \s)u ( - a) {s)dv = 0, 



S JB 



lim - / X (a) ^ a \s)^ a \s)dv = 0, (85) 

s ^°° S JB 



lim - / ^ a \s)(p {a \s)dv = 0. 



By using the Schwarz's inequality and the relation ([85]) in (|83|), we obtain 

1 r 2 

lim K c (t) - lim S c (t) = lim - / V p (a) ^ Q) (0)^ a) (2t)^. (86) 

t— >oo t— >oo t^oo 4t Jfl ^ 

When measSx ^ 0, then for u G Wi(£) , y?( a ) G Wi(fl) , the relations ((2§, QT^) and 
© imply 

/ £ u ( r\s)u^(s)dv < — / E 2W(U( S ))d« < —5(0), (87) 
Jb ^ mi Jb ^ mi 

and, by means of the Schwarz's inequality, we obtain 

l im 4 / E P {a) u^\o)uf(2t)dv} = 0. (88) 

Then, by the equations ( |8B] ) and Q8lf ) we have 

lim KM) - lim S^ft) = 0. (89) 

The relations (||) and (||) imply fl79|). 

When measY^i = 0, then, the equation (|75|) , (|76|) and (f78|) lead to 

r , 2 1,2 
1/ ; 



+y E^W+A^vmdv+y Y.p^d^dv. 

Oi = l Oi=l 



I E v i t\s)v < t\s)dv < — [ E W(V(«))di; < —5(0). (91) 

jB a=l m l J B a=1 mi 



(90) 



Since V = {v«, v( 2 \ 0«, 0( 2 )} G Wi(B) x Wi(fl) x Wi(fl) x from ([77]) 

and fl8~T|), we deduce that 
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Taking into account of the equations (P?0D and ([H]), we have 

& A / E ^ o) (0)t4 a) (2t)<fo = £ pW^^A;. 

^ o;=l ^ a=l 

With the aid of the equations (|86|) and (0), we conclude that 

lim /C c (*) = lim S c (t) + \! £ p^^dv. 
Moreover, the equations (jB2] ) and (|9~3"D lead to flSOl). 



(92) 



(93) 
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